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Abstract 

We evaluate at next-to-next-to-leading order (NNLO) the QCD corrections to the 
charm contribution from penguin diagrams to the decay K^ — )• Tr'^viy. A NNLO 
calculation is already available in the literature [1]. We provide an independent 
check of the results of non-anomalous and anomalous diagrams. We use Renormal- 
ization Group improvement and an effective theory framework to resum the large 
logarithms that appear. In the case of the non-anomalous diagrams, our results for 
the decoupling coefficients and anomalous dimensions, as well as the final numerical 
result, are in agreement with those of Ref. [Ij. In the anomalous case, analytical 
and numerical disagreements are observed. 

1 Introduction 

The rare decay mode — )■ vr+i/z/, along with — > ir^uu, plays an important role 
in flavor physics. It probes the quantum structure of flavor dynamics in the Standard 
Model (SM) [IJ-E] or its extensions [9HT3] while remaining theoretically clean. A recent 
review can be found in Ref. [H]. The cleanness of this decay is the main reason behind 
its importance, and it is due to the following: 

• It being a semi-leptonic process, the relevant hadronic operator is just a current 
operator whose matrix element can be extracted from the leading decay — t- 
7r^e~^i>, including isospin-breaking corrections |15] . 

• It is short-distance dominated. Long-distance contributions turn out to be small 
[16] . and in principle calculable by means of lattice QCD [17]. 
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Figure 1: Box- and penguin-type contributions to — t- tt+z/z/. 



As a consequence, this decay can be reliably computed with available field theoretical 
methods. Thus, the SM decay rate K'^ — )■ tt^uu, alone or together with Kl — tt^uu, 
allows for clean determinations of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. 
These decays are also very suitable for the search for new physics: a comparison of the 
determinations of various parameters from studies of CP violation in B decays with the 
determinations from the K ^ decays would provide with a clean signal of physics 
beyond the SM if substantial deviations were found. 

In the SM, the — t- ti^uu decay is mediated by box- and penguin-type diagrams, 
shown in Figure [H The leading, low-energy hamiltonian for this process can be written 
as [tlBI 



K 



Gf 



a 



eff 



y/2 2tt sin^ 6 



w 



l=e,fi,T q=t,c,u 



Here Gp, a, and sin^ 9w denote the Fermi coupling, the electromagnetic coupling, and the 
weak mixing angle, respectively. Vij are CKM matrix elements and Xq = m'^/M^. The 



contributions mediated by the exchange of virtual top quarks, encoded in the function 
X{xt), can be calculated entirely within ordinary perturbation theory. This function is 
known through next-to-leading order (NLO) P;6|, which offers sufficient accuracy. The 
contributions involving only the light u, d, and s quarks were computed in Ref. [T6] 
using Chiral Perturbation Theory, along with those of subleading, dimension-8 operators. 
When considering the charm quark contributions given by X{xc), one takes the u, d, and 
s quarks to be massless. This leads to a simplification due to the unitarity of the CKM 
matrix. Unitarity implies that 



vzvtd + v:y,, + v*y. 



'ud 







(2) 



Due to its small size, VJ* Vtd can be neglected here, and thus using Eq. (|2]) on Eq. ([T]) and 
expanding on Xc leads to 

v:y,dX{x,) + v:yudX{xu) ^ v:y^{x{x,) - x(o)) = \/:Kd(x'(o) ■ X, + . . . ) , (3) 

where higher powers in Xc are neglected. This means that in the diagrams of Figure [H and 
in their QCD corrections, one need only keep terms proportional to m^. Perturbative QCD 
effects lead to large logarithms of the form L = ln(/x^//i^), which need to be resummed 
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Figure 2: Two anomalous penguin contributions to — )■ vr+z/z/. 
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Figure 3: The one-loop top quark contribution to the coupling of the Z boson to two and 
three gluons. Diagrams obtained by the interchange of the external gluons are not shown. 

using Renormalization Group (RG) equations and an effective theory framework. This 
calculation was performed at (LO) in Refs. [H]. A study at (NLO) 

was performed in Refs. [7ll8] . which was extended to 0{a^L"'~^) (NNLO) in Ref. pj. 
The NLO calculation did not take into account the contributions from the anomalous 
penguin diagrams, such as the ones shown in Figure O These diagrams were considered 
in the NNLO calculation, but were not computed properly. A Chern-Simons operator 
was introduced, which supposedly originated from integrating out the top quark in the 
diagrams in Figure [3l However, these diagrams actually vanish in the large mt limit 
[T81IT9] . so the Wilson coefficient for the Chern-Simons operator vanishes as well, and this 
operator would make no contribution^ 

The goal of this paper is to reevaluate the penguin charm contributions to — )■ 
Ti^vv, checking the results presented in Ref. [Ij. In section [2] we briefly present our 
renormalization scheme, determined by our definition of the axial-vector current and off- 
diagonal renormalization. In section [3] we present our purely perturbative results for 
the anomalous diagrams. In section H] we perform the RG treatment of non-anomalous 
diagrams, which pervades that of the anomalous ones. Then in sections |5] and [6] we treat 
the anomalous diagrams involving a top and a bottom triangle, respectively. Finally, in 
section [7] we present and discuss our numerical results and conclude in section [HI 

^As this work was being finished, we learned that the authors of Ref. 1 had prepared an Erratum to 
their original paper in which they recalculated the contribution of the anomalous penguin diagrams. We 
discuss the Erratum in section [73] 
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2 Definition of the axial- vector current and ofF-diagonal 
renormalization 



We perform our calculations using dimensional regularization and the MS renormalizaton 
scheme. Beyond this, our scheme is determined by our choice of renormalization for the 
axial-vector current and off-diagonal corrections to the quark propagator. 

2.1 The axial- vector current 

In dimensional regularization special attention must be paid to the definition of 75. When- 
ever it appears in an open fermion line we can safely use its naive definition, but in closed 
loops we use the original definition by 't Hooft and Veltman [2D], 

(4) 

The Levi-Civita e tensor is unavoidably a four-dimensional object and thus is taken out- 
side the R-operation where a D-dimensional object can be safely considered as a four- 
dimensional one. The gamma matrices are taken as D-dimensional inside the R-operation. 
With this definition 75 no longer anticommutes with D-dimensional 7^. For this reason 
we use the definition of the axial- vector current presented in Ref. |21j, 

= (7m75 - 757/.) i^q ^J^_^^lU^U2uAqlu^lu2'lu^iJq ■ (5) 

The renormalized current is defined as = {^Z^^ Aq + Z'^Aq), where ^ = Ylq^'^ 
[2T| . The renormalization constant Z^^ cancels divergences arising from non-anomalous 
diagrams, whereas Z^ deals with anomalous ones, and thus begins at order a^. The 
coefficient ^f^ is a finite piece that is determined by requiring that the matrix elements 
of the renormalized non-singlet axial- vector and vector vertices coincide, 

(V^7m75T»^ = (V^7m^»''75 , (6) 

where is the generator of a flavor group. This effectively restores the anticommutativity 
of the 75 matrix for the non-singlet vertex, and so the standard Ward identities as well. 
This also means that the non-singlet axial-vector current defined in this way has zero 
anomalous dimension. 

2.2 OfF-diagonal renormalization 

The W and bosons can lead to flavor-changing corrections to the quark propagator, 
which in turn imply the appearance of reducible contributions to the penguin-type di- 
agrams, such as the one in Figure HI In order to avoid these diagrams, we choose to 
renormalize the left-handed doublets in the following way [22142^ . 
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Figure 4: An example of a reducible penguin-type diagram. 

where = dcos^iy + ssin^vK and s' = — d sin 9w + s cos 9w- This leads to the following 
off-diagonal counterterm in the lagrangian, 



ct 



SL U$ + gsf'Gir + -^^Y {It - sin' Ow ed) d^ + ^s^ d) 

y COS Uy/ J 

yi{Zi — Z2)sm.9w cosOw 1 (8) 

which contains the couplings to a gluon G^, with the strong coupling constant Qs and 
the S'f/(3) generator t", and the coupling to a Z boson, with the electroweak coupling 
constant Qew- I3 = (+I/27 ~l/2) and = (-1-2/3, —1/3) are the electroweak isospin and 
electric charge of up- and down-type quarks, respectively. 

We set the value of {Zi — Z2) by requiring the cancellation of off-diagonal self-energy 
corrections, 

W 4> 

y-^^^""^ + + ^ ^ . =0, 



d s d 

(9) 

plus corrections in as- This fixes also the off-diagonal gluon and Z- vertex counterterms. 

In the case of non-anomalous penguin diagrams, only the Z-vertex counterterms, like 
the ones shown in Figure [5|, survive. By definition, reducible diagrams like the one in 
Figure H] are cancelled by the insertion of the off-diagonal propagator counterterm in the 
external lines. Then we can have irreducible diagrams like the two shown in Figure O 
However, because the factor i$ in the propagator counterterm kills the denominator of one 
fermion propagator (and the factors of i introduce a sign), the sum of these two diagrams 
is zero. In the non-anomalous case we always have the same number of diagrams involving 
off-diagonal propagator or gluon counterterms, so they all cancel out against each other. 
The irreducible counterterm diagrams in the anomalous case are shown in Figure [71 Here 
the balance is broken, and one diagram with an off-diagonal gluon counterterm survives. 
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Figure 5: Off-diagonal Z-vertex counterterms at zero and one loops. 




Figure 6: Two irreducible diagrams involving off-diagonal counterterms. 




Figure 7: The counterterm diagrams for the anomalous case. 



3 Perturbative results for anomalous diagrams 

Since the anomalous diagrams where initially computed incorrectly in Ref. [1], we will 
show their calculation in detail, and begin by the perturbative results from their dia- 
grammatic calculations. We used the program qgraf [25j to generate all of the diagrams, 
and the packages q2e and exp [26||27] to express them as a series of vertices and propa- 
gators that can be read and evaluated by the FORM [2S] package MATAD 3 [25] . 

For greater convenience when resumming logarithms later on, we will express the 
results in terms of the operator Q^{fj,), defined as 

Q^(/^) = ^rT^^"7'^(l-75M' (10) 
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and write the expression for the charm contribution to HcS as 



eff 



G- - -2 
"v^27r sin^ Ow '^M^ 



9ew 



128 



where we have defined Tiu = DQy{^) and spht up the global factor comming from the 
external neutrinos. We factored out the factor 7r^/(2M^) so that the coefficient D will 
have the same normalization as the decoupling coefficients in Refs. [HE]. The relation 
between X, defined in Eq. ([1]), and D is simply X = 7r^/(2M^) ml/g'^D. 

Because of the different logarithms that appear in them, we divide the coefficient D 
in different contributions, D = D^t + -D^y,;, + D^. Here D^t stands for diagrams with a 
W boson and a top triangle, D^^i, for diagrams with W and a bottom triangle, and 
contains all anomalous diagrams with Goldstone bosons ( ^^/m^ = ) 



Dw,t = U''\li)] 12 Cf (l-4£^/,nJ(l + ^M/Mv.-^M/-c) 



(12) 



K/M^ + Kim, 



\2i 



+ 9£ 



+ 3^^/™, + 36 + 2tx^ 



(13) 



6/2 



(3 + £ 



+ 36i 



+ 75 + 2n' 



(14) 



D 



'^Mw/m.t 

+ 2 7r^ + 27 



9/2 

■^^Mw/mt 



(15) 



with / (47r) being the QCD coupling constant in the effective Uf flavor theory 

and Cp = 4/3. One can see by inspecting Eqs. f|T2l) - f|T4l) that the explicit dependence 
cancels separately in and the sum D^'^ + D^'^. 

As they are, the results from the diagrams involving the W boson are of little use 
because of the large logarithms present in them, and an RG treatment is called for. 
However, as we can see in (HM the Goldstone contribution only contains the W and top 
quark scale, which are of similar size. We won't treat this contribution when we proceed 
with the RG improvement in the effective fleld theory language. 



4 RG improvement: non- anomalous diagrams 

Part of the goal of this paper is to check the results for non-anomalous penguin diagrams 
from Ref . [l] . These diagrams are always present as subgraphs in the anomalous case. In 
particular, we will need the Wilson coefficients and anomalous dimensions related to the 
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Figure 8: Two non- anomalous penguin diagrams at NNLO. 



decoupling of the Mw scale up to NNLO, that is up to diagrams like the ones in Figure IH1 
It is useful then to present first the full RG improvement of the non-anomalous case at 
NNLO, and introduce the elements that we will need throughout. 

All penguin diagrams originate from the Z current, where the Z coupling to a fermion 
q is defined as 

Ql = VqV'^ + agA'^ = 2 (jl -2 eg sin^ Ow) V - 2 J| A'^ , (16) 
with the vector and axial-vector current 

V'^=^^ri^, and A'?=^^7'^75^, . (17) 

Since for the vector current the anomalous diagrams vanish and we want to use some 
results of the non-anomalous calculation later on, we look at the vector and axial-vector 
current separately. In the non-anomalous diagrams only the Z coupling to the four lightest 
quarks {d, u, s and c) is present, so it is convenient to define the non-anomalous current 

as 

which is the sum of all axial- vector parts of the Z current, the signs being determined by 
the prefactor — 2/|. In our renormalization scheme, the operators A^ and A'^ only enter 
in anomalous diagrams, to cancel the ones coming from A'^ and A". Since the axial- vector 
current appears only on open fermion lines in non-anomalous diagrams, we could (and 
did) use the naive definition of 75 in their computation. Using the definition in Eq. (jl]) 
leads in principle to different results: the difference between inserting one 75 or the other 
on a fermion line is of 0{e). If this line is part of a loop diagram with a pole 1/e, as in the 
penguin diagrams involving a Goldstone boson, the difference will turn finite. However, 
since this difference stems from a residue, it will be local, and thus it can be cancelled 
by adding a finite counterterm proportional to Q^. If one does this consistently at every 
order the results with a naive 75 or with 't Hooft and Veltman's 75 can be made to agree. 

We will now present the decoupling of the three heavy scales in our problem, mj , Mw, 
and rrih. Each decoupling will be assumed to take place at the appropriate scale /i close 
to the heavy scale in particular, and we will present the RG evolution that can make this 
possible afterwards. 
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(a) (b) 

Figure 9: Leading order examples of a penguin- diagram (a) and a four-quark subdiagram 
(b). 



4.1 Decoupling at the top scale 

The first scale to separate is that of the top quark, which in this case only enters through 
gluon self-energy corrections, as in Figure [HI Since we neglect all power corrections in- 
volving the top mass, this decoupling is trivial, as it is given solely by the decoupling 
relations from n/ = 6 to n/ = 5 in all the elements involved in this calculation, namely a, 
the masses, the gauge parameter ^, and the wave functions. These relations can be found 
in [30]. 



4.2 Decoupling at the W scale 

Next one needs to consider the decoupling of the W boson. Let us therefore look at the 
leading order penguin-diagram depicted in Figure [9al When considering an expansion 
of the diagram in terms of m^/M^ there are two distinct cases. The first one is when 
the loop momentum going through the u{c) line is of the order of M^y. Then the whole 
diagram shrinks to a point and can be represented by the local operator Q^, defined in 
Eq. ( ITOl) . The second case is when the loop momentum going through the u{c) line at 
the Z vertex is of the order of mc, and Mwsized loop momenta run only through a 
subdiagram containing the W line, like the one shown in Figure HH The decoupling of 
the W in the subdiagram can be described by the effective Hamiltonian 

nYd'' = ^ic^Ql + C-Q'L) , (19) 

with the definition of the effective operators (a and (3 are color indices) 
1 



± 



The combination of color indices in the definition of Q\ is chosen so that these operators 
will not mix under renormalization. The value of the coefficients will be given later on. 
After decoupling W in this subdiagram, the whole diagram is then represented by the 
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bilocal operators 



V+ 



I d'x r{ {Qlix) - Qlix)) V, V^iO) 

q=d,u,s,c 

d^x r{(Q^(a;)-Q^(x))A^^(0) 



(21) 
(22) 



for the vector and axial- vector current, respectively. The minus sign between Q'^ and 
comes from the GIM mechanism discussed in Eq. ([3]). 

Besides the operators one must also consider a set of evanescent operators 
(and the corresponding bilocal operators). We take the definitions from Ref. [1] (modified 
by a normalization factor), which we present in Appendix |Al Evanescent operators vanish 
at four dimensions, but yield non-zero contributions when inserted in loop diagrams with 
poles in e. We choose a renormalization scheme for them in which their matrix elements 
vanish and they do not mix with physical operators [3l]. The physical operators still mix 
with them, though, so the evanescents contribute to the anomalous dimensions of 
and, inserted into bilocal operators, to those of V± and A^-^. 

Putting everything together, we describe the non-anomalous penguin-diagrams with 
a vector and axial- vector current with the Hamiltonian 

H^J^ff = C+ [V^ + ] + C_ [V^ + A"!"-] + C^"^ , (23) 

Here we have omitted the normalization and neutrino factors, which are the same as the 
ones multiplying H^, in Eq. f|TT]) . The coefficients in are ( d^y^/Mw — ^^[f^w/^w] ) 



w) 



y + 2 i^i^/Mw 



-^(55 + 24C2)^/ + 26C2+ ™ 



159 -8n 



f) 



1800 



(24) 



w 



+ 4. 



■^J,w/Mw 



+ |(55 + 24C2)n/-28C2 



92443 
900 



(25) 



3 



+ 2 



34. 



-/1 14/ 



/M,^+24C2 + 33 



3 '^^lw/Mw 



24 + Uf) + il^i^^. (3408 - if^n^) 
- f V/M,. ( - 5270 - 1728 C2 + n/(161 + 72 C2) + 468 C; 



416 C4 - ^C3 + 6816 C2 - f ^/(49 + 32 (2 - 12 C; 



+ 



7995148 
675 



(26) 
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The bilocal operators V± do not mix under the RG evolution with other operators and 
their matrix elements vanish. Therefore, we don't have to consider the vector current any 
longer, its only effect being inside the decoupling constant C^^. 



4.3 Decoupling at the bottom scale 

First note that when decoupling a heavy quark, the non-singlet vector or axial-vector 
current is the same in the full and effective theory due to Ward identities p^l9|32j . As 
a consequence, the non-anomalous current is the same in the five and four flavor theory 
([A^"^]*^^) = [A^"^]*^"^)). Also, the decoupling coefficient for the Qi, operator stems purely 
from the pref actor ml/ g1, 



Qu 



(5) 



(4) 



(27) 



with ( = \n[iil/ml{iib)] 

2 



--1 



(28) 



The decoupling equations for the local operators are 



Q 



l(5) 



Ql 



l(4) 



(29) 



with 



i?_(/i,)=l + ^(a(^)(/i,)'' 



36^^,/^, + 12 + 59 



36^^,/„, + 12£^,/„, + 59 



(30) 
(31) 



These results then can be used to get the decoupling coefficients for the bilocal operators, 
whose decoupling equations read 



A 



NA 



(5) 



B4 



A 



NA 



(4) 



I rNA 

+ B±^^ 



(4) 



(32) 



with 



12- 



96- 



+ 261£^,/^, + 96C3-if^ 



B-tuif^b) 



5 («<^'(..) 



12 f 



fib/mi, 



-96- 



+ 241 96(3-^ 



(33) 
(34) 
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4.4 Decoupling at the charm scale 



At the charm scale one matches finally the operator to Qu and finds the coefficient 
X{xc), cf. Eq. ([T]). The operator multiplying X has no anomalous dimension, which 
means that X is independent of fi. In this case it is unnecessary to move to a three-ffavor 
theory, and one stays in the four-ffavor one. 



(4) 



(4) 



with the coefficients 



(35) 



C^J(/ic) =16a(^)(/i,)[l- 
C^^(/ie)=-8a(^)(/i,) 



■tMc/rric 



He /nil 



-44 



- 21 



(36) 
(37) 



After showing the decouplings at different scales, we now show how the operators can be 
evolved between these scales. 



4.5 Running to different scales 

The evolution of the operator Q^, is given by 

d 



(3^ 



where 7^ is easily determined from the renormalization constants of m^. and a^- Defining 
djji^ / /i^a = a/3(a). 



7^ = 2(7^-/3) 



3 rif 

3 \ 



-a 
9 



147 + 37 



+ 



162 



- 173259 - 18705 rif + 1535 n) + 17280 Cz^f ■ (39) 



The solution of Eq. ( |38|) is 



dz7;,(z)" 



(40) 



Considering the RG equations for the local operators Q± and A 

d 



NA 



fi'—Ql{fi) = j^Ql{fi) 



d 



' A^'^fi) = 



dii^ 



(41) 
(42) 
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one gets for the bilocal operators 



d 



The anomalous dimensions are 



7- 
and 



7+ = — 2a + a 



275267 I 52891 

1350 'V 



300 



finJ+(336 + fn,)C3 



+ a' 



12297 
50 



I 31343 ^ 260 „2 f'QQfi i 160 „ V 



7^?J 



16 a- 144 a^ + a^ 



1060082 



8 a + 208 + 



879586 
225 



225 + 144n^ + 896C3 
144n/ -64C3 



(43) 

(44) 
(45) 

(46) 
(47) 



The factor l/Qg was introduced in the definition of Q^, cf. Eq. ( ITO|) . so that the anomalous 
dimensions 7^"^ would begin at order a, and not a°. The solution of Eq. (H3l) is 



A^2^(/i) = f/±(/i,/io) A^^(/io) + f/±,.(/w,/io) g.(/io) 



where 



exp 



exp 



dz7±(z) 

a(Mo) Z /3(z) 

dz 7±(z) 



(4J 



(49) 



dz7^t(z) 



a(^to) 



/3(z) 



■ exp 



dz'7.(z')-7±(z') 



a(Ato) 



/3(z') 



(50) 



4.6 Result 



Using all the decoupling coefficients and evolution equations presented in the previous 
sections one can find the final result for the non-anomalous diagrams. We start by Eq. ( 123|) 
with Uf = 5, and evolve the operators down to a scale /if, ~ using Eqs. ( HOj) . (149|) and 
( ISUj) . We decouple the 6 quark using Eqs. ( HT^ and (1521) and continue evolving the operators 
with Uf = 4 down to a scale /i^ ~ mc, where we express everything in terms of Q^, using 
Eq. fl35|) . and then find the expression of X^^. We define the masses fhq = mq{mq) and 
the logs d-fiq/fhq = ln(yUg/mg). Further, we define the factors 



a(4)(mj 



L 



a(4)(/Xb) 



25 



M 



a(5)(/i^^) 



(51) 
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Kc comes from expressing rricifi) in terms of rric, and M and L from the evolution from 
to /ife and from fii, to /ic, respectively. With these conventions, our result reads 



X 



NA 



24 , 

J^2E J 48 r -6 



51784 r -6 /y^- 
3703 ^ 



74968 T 12 /\#12 
5819 



2888594 j^^-l , 103568 /\^-6 



103155 



18515 



149936 /\#12 
75647 



2609808 r -6 ii^-6 , 3769008 r 12 ji^l2 , T ( 317527296 nj-l 
2314375 ~'~ 3636875 ~'~ I 21490625 



212453104 /\^12 , ff _ ^496/1^-1 32 ii^-6 16 A,fl2\ 
141838125^^-' + ^Mb/nibl.-gs-^^'-' — Ta^*-^ J 



598688792 
34715625 



M- 



966244 



3772576 
40625 



486784 
21875 



M-6 + ^Mi2) +0(a2) 



(52) 



The 0{a'^) coefficient is simply too long to print here. The full result can be retrieved from 
http:/ /www-ttp.particle.uni- karlsruhe.de/Progdat a/ttpl2/ttpl2- 048/ [ We will 

give numerical results for X^^ in section [71 The expression above agrees with the result of 
the NLO calculation performed in Refs. [7f8] . correcting for the fact that in those publica- 
tions no decoupling was performed at /i^, and more importantly, that In [^c/^dfJ'c)] was 
simply changed by hand to In [/ic/mc(mc)], and so they are missing a couple terms related 
to Kc- Unfortunately, we could not readily compare the final expression of our full result 
to that of Ref. pQ, but we found agreement in the decoupling coefficients and anomalous 
dimensions leading to it (accounting in some cases for differences in normalization). 
If we re-expand Eq. (13^ around fiw we obtain 



X 



NA 



/If 2 1 4 



- ^t^w/Mw {^t^c/fhc - f ) + 1 (^Mc/m, + 27r^ + 29) 



18 



+ 



(53) 



where for brevity's sake we only show a few terms of the long, last-order coefficient. If one 
differentiates this re-expanded expression with respect to the various scales present in 
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it one finds a residual logarithmic dependence on them at the order a?, which was of course 
not present in the original, perturbative calculation. The reason is that our treatment 
is only valid up to constant terms times a^. Small logs like C-^^/md ^Mt/mj, and d^yy/Mw 
are effectively taken as constants in our treatment, and thus beyond our reach in the last 
order of Eq. (l53l) . By adding logs like these to Eq. ( 153|) one can cancel its /i^. dependence, 
which serves as a check of the correctness of our result. 



5 Top quark anomalous diagram 
5.1 Decoupling at the top scale 

The anomalous penguin diagram with a W boson and a top quark is shown in Figure llOal 
and the counterterm diagrams in Figure [3 Diagrams with the gluons attached to either 
one or both external quark lines vanish at out order (l/m°). Much like in the non- 
anomalous case, when considering the decoupling of the top quark we have two distinct 
cases. The first one is when the loop momenta of the size of mt run only through the 
subdiagram in Figure llObI The decoupling equation for this subdiagram reads 

[At] = Ct [A^] , [A^] = [A'' + A' + A' + A^ + A"] (54) 
The decoupling coefficient reads [52] ( i^t/mt = ln[/i^/mj (yUt)] ) 

Cti^t)= [a^'\^t)y[2-8i,,/m] . (55) 

The resulting penguin diagram, with the operator A^ inserted, is shown in Figure llOcI 
The second case is when all loop momenta are of the size of mt. Then the whole diagram 
llOal is treated as heavy and we get some coefficient Cu times the local operator Qi,. As it 
turns out this coefficient is zero, i.e. the subdiagram shown in Figure fTObl factorizes from 
the penguin diagram. 
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5.2 Decoupling at the W scale 

When decoupling the W we only have to treat the singlet penguin diagram in Figure llOcl 
which leads then, analogously to the non-anomalous problem, to the effective Lagrangian 

T^t.eff = Ct (C+ Al + C_ + C5 Q.) , (56) 

with the operator 

Al = -t I d'x r{{Ql{x) - Qlix)) A^{0)] (57) 

and the coefficient ( ifi^/Mw 

CUf^w) = - a^^\fiw) 8 (3 + i^^^/M^) . (58) 



5.3 Decoupling at the bottom and charm scales 



Because of the overall multiplying factor Ct, when decoupling the bottom quark we can 
neglect all terms of higher order than . This makes the decoupling equation quite trivial. 
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In the next section, when we treat the bottom anomalous diagrams, we will give this 
decoupling equation to a higher order. The singlet current in the five flavor theory contains 
the bottom quark, while in the four flavor theory only the four lightest quarks are present. 
The decoupling equation for Q^, is given in Eq. fl27|l . 
At the charm scale, the matching equation reads 



1(4) 



1(4) 
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At our order, this matching is influenced only by non-anomalous diagrams generated 
by the A'^ and A^ currents. We got rid of non-anomalous contributions from A^ and A^ 
with our off-diagonal renormalization scheme, so these operators could only enter through 
anomalous diagrams. At the charm scale, however, these diagrams would only produce 
terms beyond our precision. Nevertheless, the matching for A^. at this scale will not the 
same as the one for A^^. Since only appeared in open fermion lines we used the 
naive definition of 75 for it, but appears here through the decoupling of the current A*, 
which is inserted in a closed triangle loop. We used the definition in Eq. (|5]) for A*, and 
use it too for A? . As mentioned in section |U the different definitions of the axial- vector 
current produce different results for the penguin diagrams, which in turn leads to different 
constant parts in C\ ^ compared to the coefficients given in Eqs. ( l36l) and (ET 
( = H.^llml(^^)\ ) 
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The minus sign between these and the previous expressions comes from the different signs 
with which the currents A'^^^ enter in the definitions of and A^^. 

5.4 Running to different scales 

In contrast to A^^ the singlet operator A^ has a non-vanishing anomalous dimension 
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with flMTl 
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The solution of Eq. ( I63l) is 



dz7s(z) 



For the bilocal operator we get the following RG equation, 



A^(/xo) 



/i'^4(/i) = 7i^i(/i) + 7LQ.(/i) 



with 7I = 7± + 7^ and 



16 a- lUa^ + a^ 



7_^^ = -8 a + 80 a^ + a^ 
The solution of Eq. f l66|) is 
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with [/| and [/|^^ defined as in Eqs. ( H9|) and (|50|1 with 7± (7^^) replaced by 7I (7±^,^)- 



5.5 Result 

With the same definitions as in Eq. (I5T]) , our result for the top anomalous diagrams reads 
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Figure 11: The irreducible anomalous diagrams in the bottom case, omitting mirror 
images. 

6 Bottom quark anomalous diagrams 

The irreducible anomalous penguin diagrams with a W boson and a bottom quark are 
shown in Figure [TTl and the counterterm diagrams in Figure [71 At the order we are 
interested in there is no top quark appearing in the anomalous bottom diagrams. However, 
there are some considerations to be made before we start right away with the decoupling 
of the W boson. As mentioned in section |3l the explicit scale dependence of top and 
bottom anomalous diagrams cancels only in their sum. Since we started the evolution of 
top diagrams at the scale /it, we must choose the same starting point in the bottom case, 
lest we end up with an unbalanced /it dependence in the final result. Alternatively, given 
that d-mt/Mw ~ 1-4; we could simply perform the decoupling of the top quark at nw in the 
previous section. However, proceeding in this way effectively leads to the same problem. 
The top and W scales may be relatively large, and their difference relatively small, but 
these anomalous diagrams start already at order a^, and unlike the diagrams with the 
boson they contain also big logs involving the rric scale. Seemingly small differences can 
(and do) have noticeable effects in this case. Take the result from Eq. ( ITOl) . substituting 
/it by ixw Its /iiy dependence is, with our precision, formally the exact opposite of that of 
the result from the anomalous bottom diagrams starting at /iiy, which will be presented 
later in Eq. fl9T]) . However, numerically the dependence is about an order of magnitude 
bigger in (the modified) Eq. ( 170|) than in the bottom case, which breaks the balance 
present in the original perturbative results. Starting from /ij is formally equivalent but 
will lead to better-behaved final expressions. 

Thus, we consider the anomalous bottom diagrams starting from /it. At our order, 
the decoupling of the top quark is trivial, 

+ 0(a3) . (71) 
Next we must evolve our expressions to the scale Since the current mixes with 
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under the RG evolution, we make use of the identity 

A^ = A^ 



La^ + — a^ 



(72) 



with A^^ = A^ — l/ufA^, to express A'' through the non-mixing operators A^^ and 
A^. Our choice of renormahzation scheme for the axial current, presented in section [2?T| 
ensures that A^^ has a vanishing anomalous dimension. The evolution of A^ was shown 
in Eq. (165|) . Thus, we have that 
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The treatment of A^ is now the same as in the top quark case, and exchanging Ct by the 
coefficient from the expansion of Us, the result will be the same as in Eq. (170]) . with a 
sign coming from the different isospins of top and bottom quarks. This all leaves us with 

the operator A''{^y\r) , whose treatment we present in the following sections. 



6.1 Decoupling at the W scale 

As in the other cases we look at the different possible heavy subdiagrams to obtain the 
operators and coefficients in the effective theory. The effective Hamiltonian reads 
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and ( (!■^,^/Mw = Hf^w/^w] ) 
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We will apply the same replacement to the bilocal operator as in Eq. ( 1721) . The resulting 
operators Al^ and A\ can now be evolved separately from the W to the bottom scale. 



6.2 Decoupling at the bottom and charm scales 

For the local operators Q\. and the decoupling equations are given in the previous 
sections. For A^ the decoupling equation reads 
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with 



\n[nl/ml{^b)] ) 



2-8. 



(7J 



in complete analogy to the top decoupling in flMj) . Unlike the top quark case, however, the 
contribution from treating the whole diagram as heavy is not zero. Thus, the decoupling 
equation for the bilocal operator Aj. reads 





(5) 








(4) 




= B± Cb 


'aI 









where the coefficients B± have been defined in Eq. (!29|) . and B\. ^ read 
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At the order we are interested in, the operator A"^ for the four lightest quarks decouples 
naively ([A"^]*^^) = [A'']*^^^ + 0{a^)). Also, as in the previous section, A^ and A"^ make 
no contribution, since in our renormalization scheme they could only enter in anomalous 
diagrams that would lie beyond our precision. The decoupling of the bottom quark in a 
bilocal operator with an insertion of the current A'^ + A^ is given by 
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At our order, the difference between b'^^^^ and -B^"^, besides an obvious sign, stems from 
the fact that here we are using the 't Hooft-Veltman definition of 75 instead of the naive 
one that we used for the non-anomalous diagrams, as mentioned before. Using Eqs. fl79l) 
and (l82l) . and remembering the definition = Al' — 1/5A^ we find at last that 
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The matching at the charm scale is given in Eq. (160 
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6.3 Running to different scales 

The evolution of was given in Eq. (|69l) . A^^ has a vanishing anomalous dimension, 

(87) 



and for A^ we have 



where 
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The factor 1/5 in Eq. (!87|) comes from the definition A'"^'^ = A'' — 1/5A^. The differ- 
ence between 7^^, given in Eqs. f H6|) and P7|) . and 7^^^, comes again from the different 
definitions of 75. The solution to Eq. ( l87l) is 



(/i) = f/±(/i,/io) (/io) + ^f/^,t(/x,/xo)Q.(/io) , 
with U^^j^ defined as in Eq. (!50|) with 7^^ replaced by 7±_^. 
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6.4 Result 



Again with the same definitions as in Eq. (ISTl) . our result for the bottom anomalous 
diagrams reads 
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(5) 



the rest from 



The part proportional to [a^^\fj.t) — o-^^Hl^w)) comes from A^{ij,\y) 

1 (5) 

A'^{jjy\/) . Like in the case of the non-anomalous diagrams, the O(a^) term is too long 
;o print here, and can be retrieved from 

|http:/ /www-ttp. particle. uni-karlsruhe.de/Progdata/ttpl2/ttpl2-048/'. We give 
numerical results for X^''' in the following section. 

In the original perturbative calculation the scale dependence of the sum IS 
of order a^. The sum of the resummed results of Eqs. ( ITOi) and (19T1) satisfies this as well. 
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However, as in the non- anomalous case, the re-expanded expression has some residual 
scale dependence at order o? . It reads 



Ml 



(6) \ ^ r 2 



(92) 



The argument is the same as before. Our expansion misses constant terms times a^, and 
small logs like £^^/fh,, ^Mb/^b' ^i^w/Mw^ and i/^t/fht are effectively constants here. Adding 
terms involving these logs to Eq. (1^21) one can cancel its scale dependence and bring it to 
essentially the same form (up to pure constants) as the purely perturbative result given 
by the sum of Eqs. (fT2|) and (fT3l). 



7 Numerical Results and Discussion 



In this section we present our numerical results for the non-anomalous and anomalous 
diagrams. The analytic expression for the former is given in Eq. (1521) . and those of the 
latter in Eqs. (IT^ . flTOj) . and (191]) . We use the numerical solution of the differential 
equation relating ^^(/io) with as{^) at two loops, together with the decoupling relations 
at threshold, which we take from the Mathematica package RunDec [33j. We choose the 
following parameter values ( Th„ 



Our results are 
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Figure 12: Evolution of the relative uncertainties at each order for X^^ and X^'*"'"''. The 
bold line represents SX^j^^, the dashed one SX^^, the dot-dashed one SX^^, and the dotted 
one SX„ 



The uncertainty estimates were obtained by varying /ivk between 40 and 160 GeV, fib 
between 5 and 15 GeV and fic between 1.5 and 4.5 GeV. X^'*"*"^ represents the sum 
of top and bottom anomalous diagrams, and in this case, in addition to the uncertainty 
estimates already mentioned, we added the one fom varying fit between Tht/2 and 2mt (we 
did not touch this scale in the separate top and bottom cases as only their sum should 
be independent of it). The result for X'^ is simply the evaluation of Eq. ffl^ . times 
n'^rhl/ {2 g"^), evaluated at /i = Mw Its error, obtained by varying fj,^ as described 
above, could be reduced with an RG treatment, but its size makes that treatment rather 
pointless. 

We observe that the anomalous contributions, dominated by X^'^, are two orders 
of magnitude smaller than the non-anomalous ones, and the uncertainty of the latter 
makes them completely negligible. Indeed, the uncertainties in for X'^^ and X^'^ are of 
a remarkable size. In Figure [T2] we show the relative uncertainties of X^^ and 
at each order, that is 6X^^/\X^^\, 6{X^^ + X^to)/I^LO + X^to)\ and so on. As 
we can see, 6X^^ is completely dominated by the /Xc dependence, SXj^f^, which grows 
substantially in relative size at NLO, and keeps growing at NNLO. In the case of X^'^~^^ 
the /ic dependence still dominates and grows, although not to the same degree as in the 
non-anomalous case. 

We observe uncertainties in the pictures growing when moving from one order to the 
next. This was unexpected, as normally uncertainties should diminish when increasing 
the order of our calculation. This odd increase in fi dependence, and the dominance 
of SX^^ especially in the non-anomalous case, seem to indicate that there might be a 
problem with the perturbation series, and that it might stem from the size of ^^(/ic). At 
/X = 3GeV, we have that as — 0.252. 

We test this in two ways. First, we artificially set as{Mz) = 0.097, which leads to 
as(3GeV)~ 0.170. Second, we artificially increase the mass of the charm quark, setting it 
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Figure 13: Evolution of the relative uncertainties at each order for X^^ and X^'*"*"^, 
with an artificially small a^. The bold line represents the dashed one SX^j^^, the 

dot-dashed one 6X^^, and the dotted one 6X, 
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Figure 14: Evolution of the relative uncertainties at each order for and X^'*"*"**, 

setting nic = 3 GeV and 4.5 GeV < fic < 7.5 GeV. The bold line represents the 
dashed one SXfj_^, the dot-dashed one SX^^, and the dotted one 



to rric = 3 GeV and varying /i^ in the range 4.5 GeV< /ic < 7.5 GeV. The results are shown 
in Figures US] and HH When using the artificially smaller a^, the uncertainties in both 
X^^ and X^'*"*"* behave as we would have expected, except maybe for 6X^~^'', which 
stays stable from one order to the next. When using the heavier rric, ^Xj^^ still increases 
slightly at NLO, but its behavior is much milder than with the real m^. In the case of the 
anomalous diagrams, 6X^'^~^^ decreases and 5X^^^^ increases slightly less than before. 
From all this it seems that the low scale /ic is really straining our perturbative treatment 
of the penguin diagrams, even with the log resummation provided by the RG evolution. 
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7.1 Comparison with the previous calculation 



As mentioned in section 14.61 the final analytic expression for X^^ could not be readily 
compared with that of Ref . [1] , but we checked our agreement in the decoupling coefficients 
and anomalous dimensions that lead to it. A direct comparison between our numerical 
results for X^^ is also not possible, since only the sum of the (non-anomalous) penguin 
and box contributions (also, they chose /ic = 1.5 and fib = 5) was shown in Ref. [1]. 
Following their treatment, we computed the box diagrams and found full agreement with 
the decoupling coefficients and the anomalous dimensions they presented. However, our 
calculation was incomplete, in that we did not include the effects of the mass of the r 
lepton at NNLO. Nevertheless, choosing their values for the parameters and the scales, 
and computing what they called the "theory" error, which comes from varying the scales 
within the ranges they chose and adding the error from different treatments of (which 
they stated as ±0.001 and we did not include in our results above), we obtained for the 
quantity Pc{X) defined in their Eq. (2) the result Pc{X) = 0.376 ± 0.01, whereas in their 
Eq. (119) they presented Pc{X) = 0.375 + 0.009theory Thus, even with our incomplete 
box-diagrams calculation we have a very good agreement with their result. 

Our value for Pc{X) corresponds to X = X^ + X^^ = (9.6 ± 0.3) ■ 10"^ The relative 
error here looks much better than in all of the results in Eqs. fl94|) - fl98|) . but this is simply 
because the result is dominated by the (incomplete) box contribution, 

X^ = (9.7 ±0.2) ■ 10-\ 



X 



NA 



-(1±1) ■ w 



Our result for X^^ has changed from that of Eq. ( l94l) because of the different choice of 
lib and /ic. It is worth noting that we did not observe any of the odd behavior of scale 
uncertainties with the box diagrams. It seems that only the penguin diagrams are so 
overly sensitive. 

Unfortunately, in the anomalous case the situation is more complicated. In the Er- 
ratum of Ref. pp, the result for the anomalous diagrams is presented as APc{X) = 
1/A^(X^'* + X^'^ + X'^), where A = |K.|, for which they took the value A = 0.2248. 
As in the non-anomalous case, they chose the scales Hc = 1.5 and fib = 5. Also like in the 
non-anomalous case we could not readily compare their treatment with ours. However, 
they show an expression of their analytic result after re-expansion, and a picture with the 
evolution of APc{X) as a function of fic in the range 1 GeV < < 3 GeV at NLO and 
NNLO. The result we get for APc{X) after re-expanding Eqs. ( 1701) and ( l9Ti) and adding 
Eq. ([HD only up to NNLO reads. 
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Figure 15: APc{X) (a) and AP^(X) (b) as a function of /ic in the range 1 GeV < /ic < 3 
GeV at NLO (dashed hue) and NNLO (sohd hue). 



As mentioned before, there is a scale dependence in this expression that comes simply 
from the fact that our expansion misses terms of order times some small log, like d-^^/m^ - 
Such logs are effectively constants in our treatment and therefore missed at NNLO, and 
the fact that they are all that is needed at this order to cancel the scale dependence 
of Eq. fl99|) and bring it to essentially the same form as the perturbative result (up to 
constants) is a good check of our calculations. 
In Ref. PP the result for APc{X) is given as 

(100) 

This expression is significantly different from ours. The first thing one notices is the 
absence of logs of masses in it, which leads to the impossibility of canceling its scale 
dependence with small logs of the type ip^^jm^. One also notices the absence of /ij, which 
makes us believe that the top and W scales were considered equal. However, even if 
d^t/Mw considered small, in the perturbative result for D^'^ in Eq. (fT2|) we see that 
it is multiplied by a big log, namely imc/Mw^ ^.nd therefore one should have either fit or 
rrit in the final result. Finally, in Eq. f llOOp there is a linear if^^^/^^ that is not present 
in our case. In any case, regardless of the discrepancies of our analytic expressions, the 
size of the uncertainties makes our numerical results and theirs basically the same, both 
negligible compared to the non-anomalous case anyway. 

Regarding the evolution of APc{X) with /i^ at different orders, given that X'f' IS m- 
dependent of this (or, to our order, any) scale, we simply add X*^ as a constant both 
at NLO and NNLO. Alternatively, we define AP^{X) = AP^{X - X^). Our results are 
shown in Figure [151 Neither plot agrees with that of Ref. [1] , where the curves have the 
same shape but are translated upwards. Guided by the absence of the top scale in their 
re-expanded expression, we can find agreement with their picture if we plot AP^[X) and 
in the bottom case we neglect the running from to ^w, that is, if we compute the top 
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case starting from fit and the bottom case starting from fiw- This produces essentially 
the same results within the errors, but as mentioned before it produces an unbalanced /if 
dependence which was not present in the original perturbative diagrams. 

8 Conclusion 

In this paper we set out to check the results for penguin-type charm contributions to the 
decay — > n^uu presented in Ref. p]. It is quite an involved calculation, certainly 
deserving an independent check. We have found full agreement for all the decoupling 
coefficients and anomalous dimensions for the non-anomalous diagrams. In Ref. [1] the 
final numerical result presented was the sum of non-anomalous penguin and box-type 
diagrams, so in order to perform a comparison we also had to compute the latter. Our 
calculation missed r-mass effects, but again we confirmed the corresponding decoupling 
coefficients and anomalous dimensions presented in Ref. pQ. Incomplete as our result for 
the box diagrams was, it was enough to allow us to obtain very good agreement in the final 
numerical result, both in the central value and the estimate of theoretical uncertainty. 

We found some discrepancies with Ref. [1] in our results for the anomalous diagrams. 
Although numerically the differences get washed out by the uncertainty (and by the sheer 
smallness of the anomalous contributions versus the non-anomalous ones), the analytic 
expression of our re-expanded result differs significantly from that of Ref. p^j . 

Our results show some unstable behavior and quite a large uncertainty. The cause 
seems to be the low scale fic, which might be testing the limits of perturbation theory. 
This only applies to the penguin diagrams, however, as in our evaluation of the box 
diagrams the scale uncertainties behaved as expected, diminishing at each order and 
remaining relatively small. We were not able to find an analytic reason for the unique 
unstable behavior of the penguin diagrams. 
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A Definition of the evanescent operators 

In dimensional regularization the operators introduced in Eq. ( l20l) are not enough to 
perform the decoupling of the W boson in the sub diagram shown in Figure [9bl A set 
of evanescent operators (vanishing at d = 4) is required. In our renormalization scheme 
these operators have vanishing matrix elements, but they still contribute by determining 
the anomalous dimensions of Since we reach 0(«s) in our calculations, we need six 
evanescent operators. We take the definitions from Ref. [1] , with a different normalization 
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factor, coming from using a {V — A) current instead of left-handed fields. If we define 

Ql = {-sY{l-l,)eq){qi,{l-^,rd) (101) 
Ql = {sY{l-l,)q){qi,{l-l,)d), (102) 

where t"" is a generator of the color group, then one can define the following evanescent 
operators. 

El = {sr^'-^'a - 75)t"g) (g7,,,.M3(l - 75)t"^^) - (16 - 4e - Ae')Ql , (103) 
El = [sr^^^'-'il - 75)g) (g7m../.3(l - 75)ci) - (16 - 4e - 4e')Ql , (104) 

-(256 - 224e-^e')Q?, (105) 
El = (s7^iA«^M3M4M5(i _ (g7MiM2M3M4M5(l - 75)^^) 

-(256 -224e- i^e2)Q|, (106) 

^1 = (.V^'^^^^^^^^'^«'^^(l-75)t'^g)(g7,iM.M3M4M5Ma.r(l-75)t'^rf) 

-(4096 - 7680e)Q? , (107) 
El = (,-^m..M3M4M.M6Mr(i _ ^^)^) (g7,,,,,3M4.5M6/.7(l - 75)^) 

-(4096 - 7680e-)Q^ (108) 

Here the notation 7^i;x2...Mn = 7^17/^2 • • •7Ain has been used. The operators QI 2 are com- 
bined to produce (5± in the following way, 

Ql = Qi + ^^Q2, (109) 

71 — 1 

where ric is the number of colors. In the same way, the pairs [E^, Ef) and {E^, Eq) are 
combined to generate E± operators, which are the ones we use. The above definitions of 
the evanescent operators ensure that (5± do not mix with each other through NNLO. 
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